Condensed matter physics with trapped atomic Fermi 

gases 

H.T.C. Stoof and M. Houbiers 

Institute for Theoretical Physics, University of Utrecht, Princetonplein 5, 3584 CC Utrecht, The 
Netherlands 



Summary. — We present an overview of the various phase transitions that we 
anticipate to occur in trapped fermionic alkali gases. We also discuss the prospects 
of observing these transitions in (doubly) spin-polarized 6 Li and 40 K gases, which 
are now actively being studied by various experimental groups around the world. 

PACS 03.75.Fi - Phase coherent atomic ensemble (Bose- Einstein condensation). 
PACS 74.20.Fg - BCS theory and its development. 



1. Introduction 

After the formidable achievement of Bose-Einstein condensation in spin-polarized al- 
kali gases [1-6], the next challenge that experimentalists have already set themselves is 
to realize quantum degenerate conditions also in fermionic alkali vapors. One particu- 
lar motivation in this respect is the prediction that a gas of spin-polarized atomic Li 
becomes superfluid at densities and temperatures comparable with those at which the 
Bose-Einstein experiments are performed [7]. As a result of this experimental interest, 
the first theoretical studies of an ideal Fermi gas trapped in a harmonic external potential 
have recently appeared [8,9]. Moreover, the effects of an interatomic interaction have 
also been considered [10-12]. It is interesting to note that Oliva's calculations for atomic 
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deuterium were already performed a decade ago, even though magnetically trapped deu- 
terium had not been observed at that time. In fact, it has still not been observed, because 
the loading of the trap cannot be accomplished in the same way as for atomic hydro- 
gen [13,14]. This is presumably caused by the fact that deuterium binds more strongly to 
a superfluid helium film, that the surface recombination rate is much larger, and that the 
sample is contaminated with atomic hydrogen [15]. Fortunately, such problems do not 
arise for experiments with alkali gases and both 6 Li and 40 K have indeed been trapped 
already [16,17]. 

The most important qualitative feature of a trapped Fermi gas is that its density 
profile 'freezes' at low temperatures. This is a result of the Pauli exclusion principle and 
can be easily understood by considering an ideal gas in the trapping potential V trap (x) = 
mu) 2 x 2 /2. In general the extent of the gas cloud is determined by the classical turning 
point of the most energetic particles. If the gas is fully classical, i.e., it obeys Maxwell- 
Boltzmann statistics, these particles have an energy of a few k B T and the size of the 
cloud Rtf follows from mar R TF /2 ~ k B T, implying that R TF ^ {2k B T /muj 2 ) 1 ' 2 . We 
see that, as we lower the temperature, the size of the cloud shrinks. Moreover, if we keep 
the number of particles N fixed, the density increases. This process gradually continues 
until we reach zero temperature and the density profile becomes equal to n(x) = N5(x). 
If the gas obeys Fermi-Dirac statistics, however, the most energetic particles have in 
the degenerate (nonclassical) regime an energy that is equal to the Fermi energy ep 
and the size of the gas cloud is always given by Rtf — (2ep jmio 2 ) 1 ! 2 for temperatures 
T <C ep/k B - Comparing this also with the density profile for an ideal Bose gas, which 
in the degenerate regime consists of a large and narrow condensate peak with a width of 
about (h/mtu) 1 / 2 on top of a broad thermal background of size Rtf — (2k B T/muj 2 ) 1 / 2 , 
we conclude that the density profile of a degenerate Fermi gas is indeed 'frozen'. In 
contrast to the ideal Bose gas, the ideal Fermi gas also does not have a phase transition. 
From a point of view of condensed matter physics, an atomic Fermi gas thus appears 
much less interesting than a Bose gas. The main objective of this contribution is to argue 
that this is no longer true if there are interactions between the atoms. 

2. Interactions 

For atomic alkali gases the most important interatomic interaction is the so-called 
central interaction V c (r), which consists of a sum of the usual singlet and triplet inter- 
actions. More precisely, we have that 

(1) V c (r) = V s (r)V {s) + V T {r)V {T) , 

where and denote the projection operators on the subspace of singlet and 

triplet states, respectively. Besides this interaction that is the net result of the Coulomb 
attractions and repulsions between the electrons and nuclei of the atoms, we also have 
to consider the weak magnetic dipole-dipole interactions. Of these, the electron-electron 
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magnetic dipole interaction is most important and obeys [18] 

Here fi e is the electron magnetic moment and the tensor operator 1S obtained from 

coupling the Pauli spin matrices <Xi and cr 2 of the two valence electrons to a tensor of 
rank 2. 

It should be noted that both these interactions do not commute with the electron spin 
operators Ti(Ti/2 and therefore also do not commute with the atomic hamiltonian, which 
in a magnetic field contains both a hyperfine and a Zeeman term. As a consequence 
the central and dipolc-dipole interactions are not fully diagonal in the basis in which the 
atoms are in definite hyperfine states. This is important in principle, because it implies 
that two atoms can also collide inelastically, i.e., their hyperfine states can change during 
the collision. Together with three-body recombination events, these inelastic collisions 
in fact always seriously limit the lifetime of a trapped alkali gas. Nevertheless, we will 
in the following mostly neglect the nondiagonal parts of the interatomic interaction by 
restricting ourselves to (doubly) spin-polarized gases for which the 'good' elastic collisions 
dominate the 'bad' inelastic ones. Clearly such a restriction is a minimum requirement 
for our discussion to be also of some experimental interest. 



3. Cooper pairing 

As we will see in detail below, the most common phase transition in a weakly- 
interacting Fermi gas is due to the formation of so-called Cooper pairs. The main idea 
behind the famous Bardeen-Cooper-Schrieffer theory for this phenomenon is in fact a 
Bose-Einstcin condensation of these pairs [19]. To see how we can arrive at a mean-field 
theory for this phase transition, let us first recapitulate the mean-field (Hartree) the- 
ory for Bose-Einstein condensation in an atomic Bose gas. At zero temperature we can 
then use a variational many-body wave function in which all the atoms are in the same 
one-particle state 0(x), so 

N 

(3) *jv(xi,"-,xjv) = JJ^(xi) . 



In the language of second quantization this essentially reads 

N 



(4) \N) = 



dx </>(x)V>'(x) 



where V^( x ) is the creation operator for an atom at postion x and |0) is the 'vacuum' 
state in which there are no atoms present in the trap. Calculating the average of the 
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hamiltonian H in this variational wave function and minimizing with respect to </>(x) 
leads of course to the Gross-Pitaevskii equation [20] if N 3> 1. 

The above variational wave function has a definite number of particles. For practical 
calculations, in particular if we want to consider also nonzero temperatures or corrections 
to the mean-field theory, it is much more convenient to consider a variational wave 
function in which not the number of particles N, but instead the phase d is fixed. This 
state is given by 

(5) |0) = J2 — TTj — 10) = exp / dx 0(x)e*V(x) | 0) , 

N \J J 

as can be seen from the fact that now the matrix element (i9|^>(x)|i?) = </>(x)e*' 5 has a 
definite phase in contrast to the matrix element (N\ip(x)\N) = 0. 

It is not difficult to show that by minimizing the average energy in the state we 
again recover the Gross-Pitaevskii equation. More important for our purposes is that we 
can easily show by Fourier analysis that 

(6) lN) = j^ e -^ Wh 

which shows that the number of particles N and the phase i? are conjugate variables and 
obey the commutation relation [N, $] = i. The Heisenberg equation of motion for the 
average phase thus becomes 

< 7 » -->(§)--». 

with jU the chemical potential of the gas. We therefore recover the important Josephson 
relation 

which is also well known from the hydrodynamic formulation of the Gross-Pitaevskii 
equation [21]. Indeed, taking the gradient of this equation and using the definition of the 
superfluid velocity, i.e., v s = hV($)/m, we obtain the desired result that 

(9) 'if = 

at m 



Let us now return to Cooper-pair formation in degenerate Fermi gases. In analogy 
with Bose-Einstein condensation we can now use at zero temperature the variational 
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rfx J dx' Q , Q ,(x,x')^,(x')^(x) 



N/2 



(10) W) = ^ — l°> • 

A few remarks are in order. First, we have denoted the various hyperfine states of the 
atoms by \a). In our previous discussion of Bose-Einstein condensation we should in 
principle also have indicated the hyperfine state of the atoms and used tjj^ (x) instead of 
V> (x) in the variational wave function. However, as long as all the atoms are in the same 
hyperfine state, the precise atomic state which is trapped is unimportant from a theo- 
retical point of view and only influences the interatomic interaction, i.e., the particular 
value of the scattering length, that should be used in the Gross-Pitaevskii equation. For 
clarity we therefore suppressed the spin degrees of freedom in that case. 

Second, the fermionic creation operators anticommute. As a result the Cooper-pair 
wave function must obey <^ a>a /(x,x') = — 4> a ',a{ x ', x ), reflecting the Pauli exclusion prin- 
ciple. There are essentially two ways to fulfill this antisymmetrization requirement. If 
all the atoms are in the same hyperfine state we have a' — a and the orbital part of the 
Cooper-pair wave function is antisymmetric with respect to an exchange of the atoms. 
The relative angular momentum of the pairs must therefore be odd. In particular, we 
can have p-wave pairing. This situation occurs in doubly spin-polarized Fermi gases and 
also in liquid 3 He [22]. If we have an equal number of atoms in two hyperfine states, 
which is implicitly assumed in the above variational wave function, the spin part of 
the Cooper-pair wave function can already be antisymmetric and the relative angular 
momentum of the pairs must then be even. We now can have s-wave pairing as in ordi- 
nary superconductors. In principle, we can of course also have d-wave pairing as in the 
high-temperature superconductors, but this turns out to be extremely unlikely for dilute 
gases. 

Third, in actual applications it is again more convenient to use a variational wave 
function with a well defined phase. In this case it is obtained by multiplying the Cooper- 
pair wave function in the right-hand side of Eq. (10) with e 2 ^ and summing over all even 
values of N. We then find 

(11) |0) = cxp dx J dx' ^, Q ,(x,x')e 2 ^^,(x')^(x)^ |0) 

and exactly the same relation between the states \N) and as for a condensate of 
single atoms. The Josephson relation given in Eq. (8) is therefore also valid in this case. 
Moreover, in the state we have a nonvanishing expectation value 

(12) (i?|Va(x)Vv(x')|tf> = tfoa'toxOe 2 ** , 

which suggest that (^> a (x)^> a / (x')) is the appropriate order parameter for the Bardeen- 
Cooper-Schricffer transition, just like (^ Q (x)) is the order parameter for Bose-Einstein 
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condensation. Although this identification of the order parameter is correct, it turns out 
that it is more convenient in practice to work with the so-called BCS gap parameter 

(13) A(r; R) = Va. a ' (r)(-0 Q (R + r/2)^> a /(R — r/2)) , 

where V a ^{r) = (a, a'\V c (r) + V d (r)\a,a') is a shorthand notation for the elastic part 
of the interatomic interaction, and r = x — x' and R = (x + x')/2 are the relative and 
center-of-mass coordinates of the Cooper pair, respectively. 

Knowing the order parameter of the phase transition of interest, it is straightforward 
to obtain the corresponding mean-field theory. The main idea is first to write in the 
interaction part of the hamiltonian 

(14) = Wdx V4(x) (-^ + V trap (x) - M «) 

+ ^Y,j dX j dX ' ^'(x')^(x)^, Q '(x-x')Va(x)Va'(x') 

a, a' 

the operators ip a (x)ip a ' (x') and i/ , i(x)V , a' (x') as a sum of a mean value and fluctuations, 
and then to neglect terms that are quadratic in the fluctuations. In this manner we arrive 
at the mean-field hamiltonian 

(15) H MF = Y,j dx ^« (-^T + ~ ^ (x) ) MX) 

+ 1+ ] aa/ jdxj dx! (A(r;R)^(x')^( x )+A*(r;R)^(x)^(x')) , 

where the renormalized chemical potentials are in a good approximation (but see below) 
given by 

(16) /z' a (x) = n a - ^ J dx' V4, Q '(x-x')n Q -(x') 

and the density profile of atoms in spin state \a) obeys n a (x) = (■0^(x)-0 Q (x)). It is 
important to note that in trapped alkali gases it is indeed appropriate to have a chemical 
potential for each spin state, because the time scale for relaxation towards equilibrium 
in spin space is always much larger than the equilibration time for the spatial degrees of 
freedom. This is for example quite dramatically demonstrated by the two overlapping 
condensate experiments by Myatt et al. [4] . 

To complete the mean-field theory we should now calculate the mean values of the 
operators ip a (pc)ip a i (x') and ^ (x)i[i a i (x') in a thermal ensemble with the hamiltonian 
Hmf- This clearly makes the theory selfconsistent. To perform the calculation we write 
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the annihilation operators at the positions x and x' as 

(17) V«(R ± r/2) = fj§p M^ (TL±r/2) , 

respectively. Substituting this in the mean-field hamiltonian and neglecting gradients in 
the center-of-mass coordinate R, we find that 

(18) H MF = £ J ^ V£(k) (-^ + ^ trap (R) - /x' a (R)) ^(k) 

^ (A(k;R)^(-k)^(k) + A*(k;R)Va(k)^(-k)) . 



1 + 5 a . a > J (2tt) 

Diagonalizing the above hamiltonian by means of a Bogoliubov transformation, the spin- 
density profiles can then be calculated from 

(19) n«(R) = | ^3 (^(k)Va(k)) 

and, most importantly, the BCS gap parameter from 

/fl\c' 
— ^ t4^(k-k')(^(k')^(-k')) • 

Notice that the averages in the right-hand side of Eqs. (19) and (20) depend on the posi- 
tion R, since the mean-field hamiltonian Hmf depends parametrically on this position. 
To arrive at such a simplyfied (Thomas-Fermi) description of the inhomogencity of the 
gas we have to be able to neglect gradients of the densities n a (R) and the BCS gap 
parameter A(k;R). This is indeed true for present experiments with fermionic alkali 
gases, because the number of trapped atoms is so large that both the correlation length 
and the size of the Cooper pairs are small compared the typical length scale on which 
the trapping potential varies [11,23]. 

We refer for the details of the diagonalization to our previous work [11]. For our 
present purposes it is however important to mention that in the optimal case of equal 
spin densities, the final result of Eq. (20) is the famous BCS gap equation [19,22] 



< 2i > a * R >=-/(^ ^'< k - k '>rotr h ( 

where the so-called Bogliubov dispersion ftw(k; R) obeys 



/Mk';R) 



2k B T 



■ 2r » 2 



(22) /^(k:R.)= . «,.•:■< -|-Mk:R.)|* ■ 
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In principle, one can solve the BCS equation for any interatomic potential V a . a > (r). How- 
ever, for alkali gases the complete central interaction is usually not very well known and 
we have only information on the two-body scattering length a. We therefore would like to 
reformulate the gap equation in such a way that only this scattering length enters. This 
is achieved by noting that the gap equation is very similar to the Lippmann-Schwingcr 
equation for the two-body T(ransition) matrix. Indeed the latter reads [24] 

(23) T 2S Q ,(k,k";,) = y a>a ,(k-k")+ f -^L TW(k-k') J_^T^,(k',k";z) . 

J (27r) J z — n k' z /m 

Using the Lippman-Schwinger equation with z = 2e^(R)+zO, we can after some algebraic 
manipulation show that the gap equation is equivalent to 

< M » - -/ wfi (*Mbr» h (t) 

i 

~ft 2 k' 2 /m-2e F (R) - iO 

This result serves our purposes since the two-body T matrix is directly related to the 
two-body scattering length as we will see next. 



j A(k';R) . 



4. Doubly spin-polarized Fermi gases 

In doubly spin-polarized Fermi gases all the atoms are in the same hyperfine state. 
Although it is straightforward to generalize the following to an arbitrary state |a), we 
will in first instance restrict ourselves to the fully stretched state £} in which both the 
electron as well as the nuclear spin have a maximal projection on the magnetic field axis. 
The relevant interaction matrix element is then equal to 

(25) V u (v) = V T (r) + ^ (l - 3cos 2 (fl r )) , 

with 9 r the angle between the interatomic separation r and the magnetic field B. For 
this potential we have to calculate the two-body T matrix. Treating the weak electron- 
electron magnetic dipolc interaction in Born approximation [18] and making use of the 
fact that for the central interaction only p-waves contribute at low momenta, wc find 

(26) T c 2 f(k,k^) = ^^k.k' + ^(cos 2 (0 k _ k ,)-cos 2 (0 k+k ,)) . 

Here a is the triplet p-wave scattering length and # k ± k ' denotes the angle between the 
momentum transfer k ± k' and the magnetic field axis. 

We now have two cases to consider. Generically we expect the p-wave scattering length 
to be of the order of the range of the triplet potential and therefore &p(R)|a| to be much 
smaller than one. In that case the contribution of the triplet potential to the two-body 
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T matrix is negligible and the effective interaction between the atoms is dominated by 
the long-range dipole-dipole interaction. Due to the complicated angular dependence of 
this interaction it is not possible to solve the BCS gap equation analytically. However, 
we can nevertheless make progress by noting that the magnetic dipole-dipole interaction 
is most attractive when r is directed along the magnetic field. We thus expect that 
if Cooper pairs are formed their wavefunction ^^(r;R) is most likely proportional to 
^1,0 (r) = \/3/4:TT cos(# r ), which implies that A(k;R) = Yi i0 (k)A(R). Since this gap is 
anisotropic in the relative wave vector k, the gas is below the critical temperature an 
anisotropic superfluid, just like 3 He in the so-called A phases [22]. 

To obtain an estimate for the critical temperature we explicitly consider only the 
p-wave part of the dipole-dipole interaction, which results in the approximation 

(27) T 2 f (k, k'; z) ~ -i, W ? + ^,o)r 1 ,„ l (k)r i * m (k') . 

m 

This explicitly confirms that the dipole-dipole interaction is only attractive in the channel 
m = 0. Substituting the above two-body T matrix into the BCS gap equation and 
linearizing with respect to A(0) to obtain an equation for the critical temperature, we 
find 

(28) ldbf = J wT e F (0) V h\y2m N{k; 0) ' 

Here we used the notations V for the Cauchy principle value part of the integral and 
N(k; R) for the Fermi distribution function N(e) = (e 0e + ly 1 evaluated at ?i 2 k 2 /2m - 
£i?(R). Note that we have also used that the BCS gap equation will have a nontrivial 
solution in the center of the trap first, since the density of the gas is highest there. 
Introducing a 'scattering length' a d for the electron-electron magnetic dipolc interaction 
by means of 4:Tra d h 2 /m = — 15^oMe/2> i- e -7 fld = —^5fiofilm/8nh 2 ~ — 0.2(m/ra_y) a 
with tuh the mass of a hydrogen atom and ao its Bohr radius, the equation for the 
critical temperature becomes identical to the linearized BCS gap equation that has been 
studied previously in the context of s-wave superconductors by Sa de Melo et al. [25]. 
Using their result, we have for the critical temperature 

<2S» ^"flfiS^ 



k B TT \ 2k F {Q)\a d 

with 7 ~ 0.5772 Euler's constant. Unfortunately, the BCS transition to an anisotropic 
superfluid thus occurs at extremely low temperatures in this case and appears to be out 
of reach experimentally. For example for 6 Li at a density of n ~ 1 x 10 12 cm -3 , we have 
e F /k B ^ 600 nK and k F \a d \ ~2x 10~ 4 . 

The second case to consider appears to be more promising. As mentioned before, the 
triplet p-wave scattering length is in general too small to be able to dominate over the 
dipole-diple interaction. For 6 Li it is only —35 ao, for instance. However, we can imagine 
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that it is possible, in the same way as in the recent experiment by Inouye et al. [26], to 
(optically) trap a hypcrfine state that has a p-wave Feshbach resonance [27] and tune the 
external bias magnetic field such that the p-wave scattering length becomes large and 
negative. The two-body T matrix is then well approximated by 

(30) n B a (k,k';z) ~ l^Vk' = Mf^^Jy^jy.^) . 

m 

Substituting this into the BCS gap equation, it is not difficult to show that it is solved 
by the ansatz 

(31) A(k; R) = fcA(R) £ d* m Y hm (k) = ^A(B.)d • k , 

m 

where d m are the spherical components of a vector d that is normalized as ^2 m \d m \ 2 = 
d • d = 1. Furthermore, linearizing the BCS gap equation with respect to A(0), we find 
that the critical temperature is now determined by 

(32) ™ = / ^ k 2 7V(k;0) . 

1 ' Awa^h 2 J (2tt) 3 e F (0) - h 2 k 2 /2m 1 ' ' 

The solution can be obtained by the same methods as before and reads 

(33) Te= 8e,(0)e^/3 exp 



k B 7r "1 2(fc F (0)|a|) 3 

It is important to realize that up to this point the direction of the vector d is arbitrary. 
This reflects the rotational symmetry of the problem. However, as we have seen, the 
magnetic dipole-dipole interaction breaks this symmetry and will cause d to lie either 
parallel or perpendicular to the magnetic field, depending on precisely which hyperfinc 
state | a) is trapped. Because we are dealing with an attractive interaction we also have 
to make sure that the gas is mechanically stable. This leads to the restriction that 
n(0)|a| 3 < 5/967T, or equivalcntly (fcj? (0) | a| ) 3 < 5/576-7T 3 . Unfortunately, the latter again 
severely limits the feasibility of experimentally observing the in principle interesting 
possibility of a transition to an anisotropic superfluid. 



5. Spin-polarized Fermi gases 

We now turn our attention to the case of fermionic gases that are a mixture of two 
hyperfine states. In such a gas we can have s-wave collisions between atoms in different 
hyperfine states, which has as an advantage that interacting effects are expected to be 
much more important. Moreover, it has an additional advantage that it is now in principle 
possible to use evaporative cooling to cool the gas to low temperatures. A fermionic gas 
consisting of three hyperfine states has recently been considered as well [28] , but since it 
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does not lead to any qualitative different physics we restrict ourselves here to mixtures 
of only two spin states. We denote the two hypefme states involved from now on by | T) 
and | I). Experimentally, there are two ways to realize such a system. We can either 
magnetically trap two low-field seeking states, or use an optical trap. In the latter case it 
is presumably necessary to load the trap by precooling the gas in a magnetic trap using 
sympathetic cooling. 

5T. Positive scattering length. - If the s-wave scattering length a between two atoms 
in different spin states is positive, a first many-body effect that we have to consider is 
the phase separation of the gas into two phases with opposite 'magnetization'. Roughly 
speaking, this implies that instead of having overlapping spin densities nj (R) and n i (R) , 
the gas prefers to separate into two phases in which (almost) all the atoms are either in 
the state | f) or in the state | J.}. The driving force behind this instability is that although 
the phase separation increases the kinetic energy of the gas, this increase is more than 
compensated by the decrease in interaction energy. More precisely, the gas is stable if 
the free-energy surface F[n-\,ni\ has a positive curvature in all directions. Using that 

(34) !?>,*;*) = —, 

m 

the free energy in the degenerate regime is given by 

(35) F[n t ,nJ = J dR |(6^ 2 ) 2 / 3 ^- (nf (R) + nf (R)) + ^„ T (R)„ i( R) 

and stability of the gas requires that the spin densities obey n^(0)ni(0)a 6 < (7r/48) 2 . In 
the particular case of equal spin densities this reduces to the condition that kp(0)a < ir/2. 
Note that above we have used the T-matrix approximation to evaluate the average inter- 
action energy. Because the chemical potential fi a is equal to the derivative SF/ Sn a (R) , 
we must for consistency use the same approximation to determine the renormalized 
chemical potentials n' a (R). Therefore, we use in the following always that 

, . , , . x 4:TTah 2 

(36) ^a( R ) = Ma - 2^ - — ™«'( R ) 

m 

instead of the less accurate (only Born approximation) relation given in Eq. (16). 

Even though the interaction between the atoms is effectively repulsive, there can 
nevertheless occur a BCS pairing transition to a superfluid state due to the so-called 
Luttinger-Kohn instability [29] . Physically, the instability is a result of the fact that two 
atoms in the same hyperfine state can exchange a fluctuation (phonon) in the density 
of the other hyperfine state, leading to an effectively attractive interaction between the 
atoms involved. The p-wave transition associated with the Luttingcr-Kohn effect has 
been studied by Baranov et al. [30]. These authors obtain for the critical temperature 
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the estimate 



(37) 



e F (0) 



k B 



cxp < — 13 



( 



2fc F (0)a 



TT 



) 



2 



I 



where notably a is the s-wave scattering length. It should, however, be kept in mind 
that to be able to observe the transition we must require that the gas is mechanically 
stable and does not phase separate. As we have already seen, this implies for equal spin 
densities that kp(0)a < ir/2. Moreover, for generic values of the scattering length we 
will even have that kp(0)a < 1 at the low densities of interest. It therefore seem once 
again practically impossible to observe the above transition in a real atomic gas. 

5'2. Negative scattering length. - Our last chance of achieving a supcrfluid phase in a 
fcrmionic gas thus appears to be a spin-polarized gas with a large and negative s-wave 
scattering length. Besides the possibility of using a Feshbach resonance to tailor the 
scattering length, we can also directly make use of the anomalously large 6 Li triplet 
scattering length of —2160 a by trapping a spin-polarized 6 Li gas in a bias magnetic 
field of at least 0.05 T [31]. Substituting Eq. (34) into the BCS gap equation we see 
that the solution A(k; R) is now independent of the wave vector k and therefore equal 
to A(R). Furthermore, a linearization in A(R) gives 



Since the mechanical stability of the gas requires also for a negative scattering length 
only that fci?(0)|a| < tt/2, this expression shows that the prospects of observing a BCS 
transition in this case are indeed most favorable. 

In view of this encouraging situation we have studied in more detail the behavior of 
the spin density profile n Q (R) and the BCS gap parameter A(R) for a spin-polarized 
6 Li gas in the same magnetic trap that has been used for the Bose-Einstein condensation 
experiments with 7 Li [2,5]. The results for equal spin densities are shown in Fig. 1 and 
lead to three important conclusions. First, we see explicitly that the density profile of a 
degenerate Fermi gas is indeed completely 'frozen'. Clearly, even the BCS transition to a 
superfluid has essentially no effect on the density profile. From an experimental point of 
view this is somewhat unfortunate, because it implies that the appearance of a condensate 
of Cooper pairs cannot be observed by the same methods that have been so successful 
in the Bose-Einstein condensation experiments. We come back to this problem shortly. 
Second, for a total density of n ~ 4 x 10 12 cm~ 3 the critical temperature is about 37 
nK. In view of the achievements with bosonic alkali gases, this appears to be a density- 
temperature combination which is certainly within reach experimentally. Finally, we 
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Fig. 1. - Density distribution n-f(R) = ^i(R-) and energy gap A(R) for a 6 Li atomic gas 
consisting of 2.865 x 10° atoms in each spin state a) at T = 15 nK, b) at T = 33 nK, slightly 
below T c , and c) at T = T c = 37 nK. The left scale of each plot refers to the density and the 
right scale to the energy gap. The dotted line in c) shows the density distribution for an ideal 
Fermi gas with the same number of particles and at the same temperature. 
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have also compared the density profile of 6 Li with that of a noninteracting Fermi gas 
with the same number of particles. The strong attractive interaction between the 6 Li 
atoms evidently results in a substantial increase of the density in the center of the trap. 
Experimentally, this is a favorable effect because for a fixed number of atoms in the 
trap it enhances the critical temperature as is shown quantitatively in Fig. 2. For future 




Fig. 2. - Critical temperature as a function of the number of particles (solid line) when there are 
N/2 particles present in both spin states. The dashed line represents the critical temperature 
for a gas whose density distribution is not altered by mean-field interactions. 

convenience we mention that in a good approximation the critical temperature obeys 

,40) (£f .0.0 3 7=p{-L 36 ^ + 2.6 6 ^!} , 

with A = (2nh 2 /niksT) 1 / 2 the thermal dc Broglie wavelength and I = (fi/mw) 1 ' 2 the 
size of the harmonic oscillator ground state. It must be kept in mind that the latter 
formula can only be used for values of N x l & \a\/l that are less then about 0.74, since for 
larger values the gas is mechanically unstable and undergoes a spinodal decomposition 
first. 
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6. Discussion and conclusions 

We have argued that for condensed matter physics in trapped atomic Fermi gases, a 
spin-polarized gas with a large and negative s-wave scattering length between the atoms 
in the two different hyperfine states appears to be most promising. In view of the large 
uncertainties in the interatomic interaction potential of 40 K, the most likely candidate for 
the achievement of a gaseous BCS superfluid is at present 6 Li. However, before successful 
experiments with atomic 6 Li can be performed, some experimental problems need to be 
resolved. One problem is that in spin-polarized atomic 6 Li not only the s-wave scattering 
length but also the exchange and dipolar decay rates are anomalously large. As a result 
the gas has usually a very short lifetime. To enhance the lifetime to about 1 s at the 
densities of interest, we can either apply a bias magnetic field of about 5 T or use optical 
methods to trap two high-field seeking states. The latter solution seems to be most 
practical and is actively being persued at the moment. 

Assuming that we are able to achieve the necessary conditions for the BCS transition, 
the next problem that arises is the detection of the Cooper pair condensate. As we have 
seen, the density profile shows essentially no sign of the phase transition. Time-of- 
flight measurements, that were the 'smoking gun' for the Bosc-Einstcin condensation 
experiments, are therefore not appropriate here. Another possible signature that comes 
to mind are the frequencies of the collective modes. Because of the large s-wave scattering 
length required for relatively high critical temperatures, the collective modes are always 
in the hydrodynamical regime. They are therefore described by the (local) conservation 
laws and the Joscphson relation. In the case of equal spin densities we thus obtain the 
following set of equations. The continuity equation for a superfluid is 

(41) £~V.,. 

with the total density n = n n + n s and the total current density j = n n v n + n s v s 
consisting of a normal and superfluid contribution. In addition, Newton's law gives 

(42) ^ = - — ( Vp + n W trap ) , 

at m ' 

where p denotes the pressure in the gas. Finally, we have also the Josephson relation, 
which leads to (cf. Eq. (9)) 

(43) d -W = -~ V »- 

at m 

In principle, we must also take into account the conservation of energy. However, for 
a degenerate Fermi gas the specific heats at constant pressure and volume are almost 
equal and the continuity equation for the total energy density essentially decouples from 
the previous ones. This has important consequences, because for density fluctuations the 
Josephson relation just copies Newton's law and we must conclude that the first sound 
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modes arc not affected by the BCS transition. Of course, a measurement of second sound 
modes would be a clear signature of the transition, but this is presumably quite hard 
experimentally because second sound is primarily a temperature wave due to the fact 
that the energy fluctuations are almost decoupled. In analogy with sound attenuation in 
superconductors, it has been suggested by Fetter [32] that the damping of the first sound 
modes might be strongly influenced by the appearance of a Cooper pair condensate, but 
more work is needed to make sure whether this interesting suggestion would work. 

A property of the gas that is certainly influenced by the BCS transition is the decay of 
the gas. Qualitatively this can be easily understood from the correlator method devised 
by Kagan et al. [33] . In this approach we can relate the decay rate constant G(R; T) for 
two-body decay to the rate constant in the normal phase G(R; T c ) by 



Hence, the Cooper pair condensate enhances the decay of the gas. A quantitative estimate 
of the effect is somewhat complicated by the fact that we cannot use a pseudopotential 
to calculate the Cooper pair wavefunction ^^(R, R) from our knowledge of the gap 
parameter A(R). In any case, the increase in the two-body decay rate of the gas can 
only be used as a detection method if the gas is trapped in a magnetic trap, since in an 
optical trap the two-body decay will essentially be eliminated and the lifetime of the gas 
is determined by the rate at which photons scatter off the atoms in the gas. Because this 
is a one-atom problem, it is also not affected by the BCS transition. 

A final detection method that we would like to mention is the scattering of a beam of 
6 Li atoms from the gas cloud. Since such an experiment is quite similar to a tunneling 
experiment, it is directly sensitive to the existence of the gap parameter A(R). There- 
fore, a measurement of the angular distribution of the scattered atoms appears to be 
a promising way to get detailed information about the condensate of Cooper pairs. Of 
course, to be most sensitive we need a very cold beam. However, this is clearly not an 
impossible requirement, because the first results of such experiments with a condensate 
of 87 Rb atoms have recently been reported [34]. From a theoretical point of view, we 
are at present performing a study of the above scattering process to see how large the 
effects of the BCS transition are. It is interesting to note that it is also possible to detect 
vortices in this manner, because the atomic beam (in contrast to a laser beam) does not 
only see the core of the vortex but its complete velocity profile. This is in fact also true 
for a Bose condensate. An experiment of this sort may therefore also be of interest in 
the context of ongoing research on the properties of trapped atomic Bose gases. 
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